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Dyson-Schwinger studies of meson masses and decay constants 
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<^ ■ The masses and decay constants of the light mesons are studied within a ladder-rainbow 
■ truncation of the set of Dyson-Schwinger equations using a model 2-point gluon function. 
The one phenomenological parameter and two current quark masses are fitted to reproduce 
f n , m n and mx- Our results for fx, and for the vector mesons p, 4>, and K* are in good 
| agreement with the experimental values. 

O ■ 1- DYSON-SCHWINGER EQUATIONS 

00 , 

The Dyson-Schwinger equations [DSEs] form an excellent tool to study nonperturbative 
! aspects of quark propagators and the formation of bound states |TJ. The approach is 
consistent with quark and gluon confinement f2[ and generates dynamical chiral symmetry 
breaking. It is straightforward to implement the correct one-loop renormalization group 
behavior of QCD ||, and obtain agreement with perturbation theory in the perturbative 
region. Phenomenological models based on the DSEs {§] have recently allowed significant 
progress not only for light mesons, but also for heavy mesons |J. 

The DSE for the renormalized dressed-quark propagator in Euclidean formulation is 

Sip)" 1 = Z 2 i 1 -p + Z 4 m q ( f i) + Z 1 J A ^g 2 D^(p-q)^S(q)r:(q,p), (1) 

where D^ u {k) is the renormalized dressed-gluon propagator, T®(q;p) the renormalized 
dressed- quark-gluon vertex, and J A denotes a translationally-invariant regularization of 
the integral with A the regularization mass-scale. The general solution of Eq. (|l|) has the 
form Sip)' 1 = 27 • pA(p 2 ) + B(p 2 ), renormalized according to 
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sip)- 1 



iy ■ p + m q (p,) , (2) 



with m q (fi) the current quark mass at the renormalization scale fi. The renormalization 
constants Zi, Z 2 and Z4 generally depend on fi and A, but not on the quark flavor. 

The renormalized, homogeneous Bethe-Salpeter equation [BSE] for a bound state of a 
quark of flavor a and an antiquark of flavor b having total momentum P is given by 

Tt(p; P) = J A ^K(p, q- P)S a (q + vP^tiq; P)S b (q - fjP) , (3) 

where 77 + fj — 1 describes momentum sharing, (p; P) is the Bethe-Salpeter amplitude 
[BSA], and M specifies the meson type: pseudoscalar, vector, axial- vector, or scalar. The 
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kernel K operates in the direct product space of color and Dirac spin for the quark and 
antiquark and is the renormalized, amputated qq scattering kernel that is irreducible with 
respect to a pair of qq lines. This equation defines an eigenvalue problem with physical 
solutions at the mass-shell points P 2 = —m 2 with m being the bound state mass. 
We use a ladder truncation for the BSE 

Ki:(p,q^P)^-g((p-q) 2 )D^( P -q)i [ j 1 , ] j ® f^J , (4) 

where D f ™ c (k) is the bare gluon propagator in Landau gauge. Together with a rainbow 
truncation — > ■y u \ a /2 for the quark DSE, Eq. (||), this preserves the axial-vector 
Ward-Takahashi identity ||, which ensures that in the chiral limit the ground state 
pseudoscalar mesons are massless even though the quark mass functions are strongly 
enhanced in the infrared. Thus pions are understood as both the Goldstone bosons 
associated with dynamical chiral symmetry breaking and as qq bound states [0,0. 
For the "effective coupling" G(k 2 ) we employ the Ansatz || 



1 + k 2 /A 2 QCD 



Hk 2 ) , (5) 



with F{k 2 ) = [1 - exp{-k 2 /[Am 2 ])}/k 2 , r = e 2 - 1, and 7m = 12/(33 - 2^/). This 
Ansatz preserves the one-loop renormalization group behavior of QCD; in particular, 
one obtains the correct one-loop QCD anomalous dimension of the quark mass function 
M(p 2 ) := B(p 2 )/A(p 2 ) for both the chiral limit {m(fi) = 0) and the explicit chirally bro- 
ken case {m(n) ^ 0). The first term of Eq. (^|) implements a strong infrared enhancement 
in the region k 2 < 1 GeV 2 which is phenomenologically required for sufficient dynami- 
cal chiral symmetry breaking and for confinement 0. We use m t = 0.5 GeV, Nf = 4, 
^qce> 4 = 0.234 GeV, lo = 0.4 GeV, and a renormalization point /i = 19 GeV. 

2. MESON BETHE-SALPETER AMPLITUDES 

The BSA Tps(q; P) of a pseudoscalar meson has the general form 

T PS (q;P) = l 5 (tE + 1 -PF + 1 -qG + q^P v H) , (6) 

with the invariant amplitudes E, F, G and H being Lorentz scalar functions of q 2 and 
q ■ P; P 2 = —m 2 PS is the fixed, on-shell, meson momentum. For unflavored mesons that 
are charge conjugation eigenstates, such as the ir, there is an additional constraint on the 
BSA to obtain a specified C-parityQ Of the four covariants given in Eq. (^), G is odd 
under C, the others are even under C. The only remaining quantity that can produce a 
desired uniform C-parity is q-P which is odd under C. Thus a pion will have an amplitude 
G that is odd in q ■ P while the other amplitudes are even in q ■ P. Flavored mesons like 
kaons, will have amplitudes that are neither even nor odd in q ■ P. 



x We do not discriminate between up and down quarks, nor incorporate electromagnetic corrections; 
therefore the BSA for 7r ± is equal to that for 7r°. 
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Massive vector mesons are transverse, and the general form of the BSA T^(q;P) can 
be expressed in terms of eight independent Lorentz covariants T^(q; P) with coefficients 
Fi(q 2 , q ■ P; P 2 )- The choice for the covariants T l {q; P) to be used as a basis is constrained 
by the required properties under Lorentz and parity transformations, but is not unique; 
a convenient orthogonal basis is given in Ref. |§. For the other mesons, scalar and axial- 
vector, covariants can be defined similarly. However, we expect the ladder truncation 
used here to be reliable for pseudoscalar and vector, but not for scalar mesons |||K|- 

The electromagnetic decay mediated by a photon (e.g. p°, lu, (ft), and the leptonic decay 
mediated by a W-boson (e.g. p ± , K ), are described by the vector decay constant P,|8|JTT[| 



f v m v e£\P) = (0\q%q a \V ab (P,\)), (7) 
where ejp is the polarization vector of the vector meson. This is analogous to 



fpsP, = (0\q% l5 q a \PS a \P)) , 

= Z 2 N C / A 704Tr D [ W S a (q + r)P)T%(q;P)S b (q-fjP)], (8) 

for the pseudoscalar decay constant corresponding to f w = 131 MeV under the normaliza- 
tion convention used here. The vector decay constant can be expressed as a loop integral 

fvmv = ~y^J J^Tr D [ r ,S a (q + V P)rf(q- 1 P)S b (q-fjP)}. (9) 
Both Eq. (H) and(||) are exact if the dressed quark propagators and the BSA are exact ||[7||- 

3. RESULTS FOR PSEUDOSCALAR AND VECTOR MESONS 

After insertion of the general form of Tps(q; P) or T^(q; P) in terms of a covariant basis, 
the homogeneous BSE Eq. (^), can be reduced to a set of coupled integral equations for 
the scalar functions F? b {q 2 , q ■ P\ P 2 ) 

F?(p 2 ,p.P;P 2 ) = J A -^K^(p 2 ,p-P-,q 2 ,q-P;p-q-,P 2 )Ff(q 2 ,q-P;P 2 )(10) 

This equation is solved numerically without further approximations as an eigenvalue prob- 
lem to give the mesons masses; the decay constants follow from Eqs. (|8p and (0). 

The parameter D and the current quark mass m u / d (fi) are fixed by fitting and f n . 
Next, the strange quark mass m s {fi) is determined by a fit to tuk- With D = 0.93 GeV 2 
and quark masses m u /d{^) = 5.5 MeV, and m s (fi) = 130 MeV scaled to /i = 1 GeV, we 
obtain a good description of the n, p, K, K* and <fi masses and decay constants, see 
Table |1|. Our results are not very sensitive to the details of the gluon Ansatz M , similar 
results are found with u = 0.3 GeV and uj = 0.5 GeV, as long as D and the quark masses 
are refitted. Also r n is in excellent agreement with the experimental value ]T^| . 



For small current quark masses, the pseudoscalar meson mass grows like m 2 PS oc m g (p), 
in agreement with the Gell-Mann-Oakes-Renner relation. For larger quark masses, 
m q (p) > 3m s (/i), the meson mass grows linearly with the current quark mass fl3[ 



as 



expected for heavy mesons. The decay constants also increase with increasing quark 
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Table 1 

The pseudoscalar and vector meson masses tum and decay constants ju in GeV. 





m (calc.) 


/ (calc.) 


m (exp.) 


/ (exp.) 


pion 


0.138 


0.131 


0.1385 


0.1307 


kaon 


0.497 


0.155 


0.496 


0.160 


P 


0.742 


0.207 


0.770 


0.216 


K* 


0.936 


0.241 


0.892 


0.225 





1.072 


0.259 


1.020 


0.237 



masses, as can be seen from Table |l|. However, for heavy mesons it follows from the 
general form of the BSA, together with the normalization condition, that both the vector 
and the pseudoscalar decay constant decrease with the quark mass as /m oc 1 / ^m M ||. 

For mesons that are charge conjugation eigenstates, the dependence of the amplitudes 
Fi(q 2 , q ■ P; P 2 ) on q ■ P is minimal. For flavored mesons there is a much stronger angular 
dependence. Our results for physical observables, such as the mass and decay constant, 
are independent of the momentum sharing and of the regularization scale A, as long as all 
relevant covariants and the full angular dependence are included in the calculation. 
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